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We present perturbative analytical results of the Landau level quasienergy spectrum, autocorrela¬ 
tion function and out of plane pseudospin polarization for a single graphene sheet subject to intense 
circularly polarized terahertz radiation. For the quasienergy spectrum, we find a striking non trivial 
level-dependent dynamically induced gap structure. This photoinduced modulation of the energy 
band structure gives rise to shifts of the revival times in the autocorrelation function and it also 
leads to modulation of the oscillations in the dynamical evolution of the out of plane pseudospin 
polarization, which measures the angular momentum transfer between light and graphene electrons. 

For a coherent state, chosen as an initial pseudospin configuration, the dynamics induces additional 
quantum revivals of the wave function that manifest as shifts of the maxima and minima of the 
autocorrelation function, with additional partial revivals and beating patterns. These additional 
maxima and beating patterns stem from the effective dynamical coupling of the static eigenstates. 

We discuss the possible experimental detection schemes of our theoretical results and their relevance 
in new practical implementation of radiation fields in graphene physics. 

PACS numbers: 78.67.Wj, 71.70.Di, 72.80.Vp 


I. INTRODUCTION 

The dynamical control of the transport properties of 
Dirac fermions in the condensed matter realm is currently 
an intense research topic. These Dirac fermions have 
been shown to emerge as the low energy excitations of 
two-dimensional systems with a honeycomb lattice struc¬ 
ture as it occurs in grapheneir— . Recent theoretical^— 
and experimental^ works have discussed the role of radi¬ 
ation fields in the manipulation of the transport proper¬ 
ties of monolayer graphene samples. By focusing on the 
Terahertz frequency regime particular attention is paid 
to the tunability of the induced band gaps. In addition, 
the possibility of generating topological insulating behav¬ 
ior was theoretically put forward both in the static^ and 
dynamical regimes— 7— 

In presence of a perpendicular quantizing magnetic 
field B = Bz, the static spectrum of single layer graphene 
posseses a \J~B field dependence which strikingly con¬ 
trasts the linear B dependence for conventional non rela¬ 
tivistic 2DEG—i. In addition, the n = 0 Landau level (LL) 
has only one sublattice occupied at each Dirac point. 
Considering the LL scenario and a topological contribu¬ 
tion given by an excitonic gap the authors of reference^ 
predict the appearance of Rabi oscillations when the sys¬ 
tem’s initial quantum state is prepared by means of a 
short electric pulse and the subsequent dynamics is con¬ 
trolled by the oscillations between the dynamically cou¬ 
pled LL. 

In this work we theoretically analyze the dynamical 
manipulation of the LL structure of charge carriers on 
suspended monolayer graphene when a periodically driv¬ 
ing radiation field is applied perpendicular to the sam¬ 


ple. A similar set up was proposed in reference^, where a 
Gaussian laser pulse is introduced via the dipolar approx¬ 
imation. In that work, the authors discuss the dynamics 
of Zitterbewegung which is described in terms of the radi¬ 
ation field emmited by the accelerated charge carriers in 
graphene. In our proposal we consider a continuous laser 
field and thus make use of Floquet’s theorem to recast the 
dynamics in an explicitly time-independent fashion but 
without need to resort to infinite-dimensional Fourier¬ 
mode expansion. Our approach has the advantage of 
providing an analytical description of the driven evolu¬ 
tion of relevant physical quantities as the pseudospin po¬ 
larization which is a measure of the angular momentum 
exchange among the charge carriers and the radiation 
field^. In the following we explicitly show how our semi 
analytical results are relevant at high values of the quan¬ 
tizing magnetic field B , where the coupling to the radi¬ 
ation field leads to non trivial qualitative modifications 
of the dynamical behavior of relevant physical quanti¬ 
ties within the perturbative regime. To the best of our 
knowledge, our work constitutes the first approach to the 
photoinduced modulation of the Landau level structure 
in single-layer graphene in presence of an intense and con¬ 
tinuously applied laser field. However, we would like to 
mention that a recent worbi^ has addressed the role of 
intense radiation field. Yet, the authors of this work con¬ 
sider a quantized radiation field and do not address the 
Landau level structure scenario in single-layer graphene. 

The paper is organized as follows. In section II we present 
the model and summarize the perturbative results for the 
quasienergy spectrum. In section III we study the dy¬ 
namics of the autocorrelation function and pseudospin 
polarization for two initially prepared states. First we 
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consider an eigenstate of the static Hamiltonian that has 
vanishing pseudospin polarization in the static regime. 
Next, we present the same analysis for an initially pre¬ 
pared coherent state, with static finite pseudospin polar¬ 
ization and highlight the main differences with respect to 
the other initial configuration. In section IV we discuss 
our main results and in section V we present conclud¬ 
ing remarks and argue on the experimental implemen¬ 
tation of our proposed theoretical setup. Finally, in the 
appendix we summarize some mathematical calculations 
arising during the perturbative analysis. 

II. MODEL 

In this section we focus on the low energy properties 
of non interacting spinless charge carriers in a suspended 
monolayer graphene subject to a perpendicular, uniform 
and constant magnetic field B = Bz. The dynamics is 
governed by Dirac’s Hamiltonian. In coordinate repre¬ 
sentation it reads 

H v (r) = v F (r]ir x a x + TT v a y ), ( 1 ) 

where v F ~ 10 6 m/s is the Fermi velocity in graphene. In 
addition, the canonical momenta nij = pj + eAj ( j = x, y ) 
contain the vector potential (V x A = B), — e is the 
electronic charge (e > 0), and y = ±1 describes the 
valley degree of freedom. Using the definition of the 
magnetic length l jj 2 = eB/K and the cyclotron energy 
tvjj c = \[2 v f H/Ib the Hamiltonian Eq. (JTJ at each K (K’) 
Dirac point, which corresponds to 77 = +1 (77 = — 1), can 
be written in the form 


for n/0. The zero energy eigenstate (n = 0) is given in 
each case by 

It’D = (10)) < 7 > 

te 1 } = ( 'o ) ■ ( 8 ) 

Due to time-reversal symmetry, we have TH + \T = H—\. 
Let us now consider the effect of intense circularly po¬ 
larized terahertz electromagnetic radiation, incident per¬ 
pendicularly to the sample. We assume that the beam 
radiation spot is large enough compared to the lattice 
spacing so we can neglect any spatial variation. Accord¬ 
ing to the standard light-matter interaction formulation, 
the dynamical effects of a monochromatic radiation field 
incident perpendicular to the sample can be described by 
means of a time-dependent vector potential 

£ 

A (t) = — (cos ut, 5 sin u/t) , (9) 

u> 

where £ and uj are respectively the amplitude and fre¬ 
quency of the electric field given in turn by the stan¬ 
dard relation £{t) = —d t A(t). In addition, we are us¬ 
ing S = +1 (5 = — 1 ) for right (left) circular polar¬ 
ization. We are using circular polarization because it 
has been shown to provide the maximal photoinduced 
bandgap^. Starting from the ordinary dipolar interac¬ 
tion term —ep ■ A (t), introduced to the Tight-Binding 
Hamiltonian via the Peierls substitution, we can evalu¬ 
ate the effects of the driving at each Dirac point as 

V v = ev F [rja x A x (t) + Sa y A y (t)}, (10) 
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where the annihilation and creation operators are defined 
by standard relations as 


a = 




y j anda t = (*•*+»%) 


V2 V2 

The eigenenergies of the Hamiltonian 0 are then 


E% v = spy/ntujjc (4) 


with s = ±1. Positive (negative) values of sy represents 
the conduction (valence) band at each Dirac point. In ad¬ 
dition, the integer quantum number n = 0 , 1 , 2 ... labels 
the Landau level (LL) structure of monolayer graphene. 
Using the eigenstates | n) of the operator a^a, the corre¬ 
sponding eigenstates | <Pn V ) read 


} ~V 2 \ l n > ) 
i^.-u -_L( - s l n > ^ 


( 5 ) 


which explicitly reads 

V v = £r](ux cos cut + rjSa y sin ut). ( 11 ) 

with the effective coupling constant £ = ev F £/ui. This 
makes the total Hamiltonian 

H r ,(t) = H v + V v (t), (12) 

periodic in time H v (t + T) = H v (t ), with T = 2ir/u> the 
period of oscillation of the driving field. Therefore, if 
we focus on the K Dirac point (77 = 1), the physics at 
the K’ Dirac point (77 = —1) can be easily found by the 
substitutions £ —> — £ and w —> —u). 

Thus, let us focus on the K point physics and afterwards, 
we can make the necessary substitutions. In order to sim¬ 
plify the notation we set H + 1 = H 0 and V+i(f) = V(t). 
Hence, defining rising <r+ and lowering er_ pseudospin 
operators by the standard formulas 

(J x T iCTy 


the time-dependent interaction potential can be rewrit¬ 
ten as 


( 6 ) 


V(t) = £(e- i 5 “V+ + e iSut a-), 


(13) 





3 



FIG. 1. (Color online) Approximate mean energies for the driven scenario (continuous lines) as function of the quantizing 
magnetic field B. The dotted lines represent the undriven spectra for the corresponding LL. We notice a level-dependent energy 
gapped that leads to non trivial behavior of physical quantities, as discussed below (see main text). 


Now we invoke Floquet’s theorem which states that 
the time evolution operator of the system induced by a 
periodic Hamiltonian can be written in the formic 

U{t) = P{t)e ~ lHFt / h , (14) 

with P(t) a periodic unitary matrix and Hp a time- 
independent dynamical generator referred to as the 
Floquet Hamiltonian. The eigenvalues of the Floquet 
Hamiltonian Hp represent the quasienergy spectrum of 
the periodically driven system. Typically, in order to 
solve for the quasienergy spectrum, one can expand each 
term of the time-dependent Schrodinger equation in 
Fourier space and numerically solve an infinite eigenvalue 
problem. Instead, we will take a perturbative approach 
as discussed below. 

Accordingly, for our problem we can find approximate 
solutions to the dynamics by modifying slightly the ana¬ 
lytical strategy presented ii*i£. Then, one finds that the 
excitation number operator N a , defined as 

= (a t a+^)l + y, (15) 


which commutes with the Hamiltonian Hq and satisfies 
the eigenvalue equation 

N aWn) = n\V S n)- (16) 

N a generates a time-dependent unitary transformation 
|’F(t)} = P(t)|$(£)) given as 

P{t) = exp(— iN a 5u)t), (17) 

such that the time-dependent Schrodinger equation 

ihd t \y(t)) = H(t)\*(t)) (18) 

can be transformed with a time-independent operator Hp 
governing the dynamics of the problem 

mm)) = H F \$(t)), ( 19 ) 

where Hp and |$(t)) are the Floquet Hamiltonian and 
Floquet eigenstate, respectively. Doing the explicit cal¬ 
culation, Hp is found to be given by 

Hp = H 0 - N a 8huj + £a x . 


N c 


( 20 ) 
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In the following we focus on recent experiments in the 
far infrared frequency domain^ for which hio ss 10 meV 
and we consider values of the electric field intensities 
£ ~ 0.15 MV/m. Then one gets for the coupling con¬ 
stant £ « 10 meV which, for frequencies oj in the tera¬ 
hertz domain leads to £ ~ tko. This value is an order 
of magnitude smaller than the Landau level separation 
hu c ~ 116 meV, for B = 10 T. For larger frequencies 
and stronger magnetic field intensities, the ratio £,/tujj c 
tends to be smaller. Therefore, we can perform a per¬ 
turbative treatment in the effective coupling parameter 
A = £/hu) c < 1. 

We should remark that although our radiation field is 
intense, it is one order of magnitude smaller than the nu¬ 
merical estimates used in reference^ for which one gets 
£ ~ 1.5 MV/m. Yet, the experimental setup used in 
reference^ consisted of infrared radiation field with power 
P = 20 mW, sample areas equal to A = 3 x 3mm 2 and 
d = 5x 5mm 2 ; thus, one gets electric field intensities of 
order £ ~ 1 kV/m, which in turn leads to £ ~ 3 meV. 
In this manner, our perturbative results would allow for 
an analytical treatment of future experimental extensions 
of the work described in£, in case they would include a 
quantizing magnetic field in their study. It would also al¬ 
low for larger values of the radiation field intensity with 
novel photoinduced features in the Landau level struc¬ 
ture of single- layer graphene as it is described in the 
following. 

For this purpose we transform the Hamiltonian (1201) as 
H = , where we have introduced the 

antihermitian operator /_ = a'cr_ — ao+. Evaluating up 
to first order we get 

H*H f + ±\I-,H f ]. (21) 


Evaluation of the commutator gives (in the appendix we 
summarize the explicit derivations) 


\I ~, Hp] — — buj c 2!V 0 T A(o/ -t- cl) 


( 22 ) 


and defining the shifted operator b = a + A one gets, to 
leading order in A, the effective Hamiltonian 


H = huj c (b^a- + b(j + ^J — ShuNb — £,Nb<J z , (23) 


where we have introduced the shifted number operator 

N b = tfb + ?-^. 

In equation (1231) we have neglected the additive higher 
order terms 


AV = hco A(b f + b) - A 2 hu], (24) 

which can be dealt with by higher order perturbation 
theory. Thus, the approximate quasienergies are found 
to be given as 

e s m = sHu} c y/m \/1 + mA 2 , mod Hlo, (25) 


which can be rewritten as 


e s m = s\Jm(fojj c ) 2 + (?n£) 2 , mod hw. (26) 


Then, to this order of approximation, all quasienergies 
corresponding to the m ^ 0 LL are shifted, whereas 
the m = 0 remains insensitive to the radiation field. 
This shift of all but the m = 0 quasienergy LL spec¬ 
trum agrees with the result reported in reference^ for 
bilayer graphene. Yet, a more detailed derivation by 
means of second order perturbation theory shows that 
there is a small 0(A 4 ) energy correction due to first non¬ 
diagonal terms in equation (l24l) which couple all adjacent 
LL. This higher order corrections could be important at 
low quantizing magnetic fields for which the condition 
A = £/hu} c ~ 1 could be satisfied. 

Let us now analyze some physical consequences of the ra¬ 
diation field on the Landau level structure of monolayer 
graphene with focus on the interplay among the quantiz¬ 
ing magnetic field and the light-matter interaction. To 
begin with, we notice that the quasienergies are defined 
up to multiples of w; therefore, a better physical charac¬ 
terization of the energy spectrum for the driven system 
is provided by the mean energies^ 


de s 


^ rm ^ 'YY1 ^ 


duj ’ 


(27) 


which are invariant under e s m —> e s m + IHui, for l being an 
integer. Doing the explicit calculation the mean energies 
are found to be given by the expression 


= s 


m{tvuj c ) 2 + 2m 2 £ 2 ^ 
\/m(Hu c ) 2 +?n 2 £ 2 ) 1 


(28) 


where we remember the definition of the effective 
coupling to the radiation field as £ = ev F £/oj. 

As can be seen in FIG. 1, these mean energies are 
plotted as function of the quantizing magnetic field H, 
for different values of the Landau level index changing 
the effective coupling £. We notice that, at intermediate 
light-coupling strength, the energy resolution of these 
levels becomes much better and could experimentally 
be tested for not so large quantizing magnetic fields B. 
Moreover, we find that to this order of approximation 
the LL become gapped, with the striking feature that 
the photoinduced gap is level dependent. These gap 
openings appear except for the m = 0 level which, as 
discussed before, remains insensitive to the radiation 
field. 


Below we will deal with the photoinduced dynamical 
features and therefore, we give the corresponding nor¬ 
malized Floquet eigenstates for m / 0 

| 0=(' J AW 1> ) > « 

where we have defined the coefficients 


r 

J m 


smt; 


2 e n 


(30) 
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with e m = |e£j. In addition, the zero energy eigenstate 
(to = 0 ) is still given by |^>o}> i- e., 

I^o) = ( |J) ) • (31) 


For the chosen initial state we find cr z (t,0) = 0. When 
expression © is evaluated we find after some algebraic 
manipulations 


Vz{t,Q 



sin 2 e m t/h. 


(37) 


III. PSEUDOSPIN AND AUTOCORRELATION 
FUNCTION DYNAMICS 

Now that we have found the approximate Floquet 
eigenstates and quasienergies, we explore other dynam¬ 
ical features of the driven LL configuration by evaluat¬ 
ing the mean values of the pseudospin polarization op¬ 
erator and its relation to the autocorrelation function 
dynamics^. 

For this purpose, let us first assume that the system is 
initially prepared in an eigenstate of the Hamiltonian Hq 

l*(0)> = IO, (32) 

with m ^ 0. In the Floquet basis (l29l) . the initial state 
is written as 

kj= E D ™\ (33) 

s'—=hs 


We note that for vanishing values of the coupling to the 
radiation field £ —;> 0 one has a z (t, 0) = 0. Therefore, 
once the electromagnetic field is present the pseudospin 
oscillations are a manifestation of the angular momen¬ 
tum exchange among the radiation field and the charge 
carriers in graphene^. As discussed in the case of the 
quasienergies, we could better quantify the effects of the 
driving field by evaluating the average 

K) = ? L (38) 

with the period of the radiation field given as T = 27r/w. 
Then we get the expression 

(er z ) = (i _ sine 2e m T/K) , (39) 

^ m 

with sine a; = sinx/x. 


where the expansion coefficients are given by 

D+ ss 'fm') ■ (34) 

Taking into account that the | t/A) states are degenerate 
eigenstates of Nf, , the unitary operator e ~ lN b^ t would 
just contribute a phase e~ lujmt . Using the Hamiltonian 
I© the evolved state can be written as 

!*(*)>= E D^{t)e~ is '^ t/h \^i), (35) 

S'—±S 

where we have introduced the time-dependent coef¬ 
ficients D™' (t) = Let us then consider 

the dynamics for the autocorrelation function and 
pseudospin polarization cr z (t) operators from which we 
can respectively infer the feasibility of manipulating 
the polarization state of the sample and the dipole 
moment radiation emmited by the driven sample. It 
is important to remark at this point that studying 
the out of plane pseudospin polarization is a means of 
detecting the angular momentum exchange between the 
Dirac fermions in graphene and the circularly polarized 
radiation field, as can be inferred from the discussion 
in the recent literature about the role of er- in describ¬ 
ing the total angular momentum content of the system- 14 . 

With these ideas in mind, we first begin by evaluating 
the pseudospin polarization 

(36) 


In figure FIG. [2] we plot the behavior of (a z ) as a func¬ 
tion of the coupling strength £. At small £ = 0.25 Hco we 
see that the closer the state is to the to = 0 LL the role 
of the radiation field in modifying its pseudospin polar¬ 
ization is less relevant and this is correlated to the fact 
that, within this regime, the gap openings seen in the 
mean quasienergy spectrum are not so noticeable at dif¬ 
ferent values of the to LL index. In addition, at large 
values of the quantizing field, B = 10T, the value of 
the pseudospin polarization is almost the same at each 
value of the to LL index. However, we see in the sec¬ 
ond upper panel that already at intermediate values of 
the relative coupling strentgh £ = 0.5fiw that at large 
values of the quantizing field B = 10T one can discern 
among the different LL polarization value which serves to 
separate each level contribution to this quantity dynam¬ 
ical behavior. In the lower two panels we show that this 
physical picture becomes more obvious at larger values 
of £, with perfect separation for the contribution from 
the to = 1 and m = 2, levels. Yet, the corresponding 
to = 3 and to = 4 contributions to the pseudospin polar¬ 
ization remain degenerate at large B values. One would 
expect that this accidental crossings seen as a degeneracy 
would be removed by including higher order perturbative 
contributions, since then the coupling among nearby lev¬ 
els would lead to additional splittings in the quasienergy 
spectrum that would in turn also lift this accidental pseu¬ 
dospin degeneracy. 

In order to discuss more general dynamical features, 
one should consider a superposition state. However, in 
order to do so we must take into account the feasibility 


Vz(t,0 = <^(t)M^(t)). 
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FIG. 2. (Color online) Approximate mean pseudospin polarization (<r z ) for the driven scenario as given by equation (37), 
plotted as a function of the quantizing magnetic field B. We have selected the first four LL m = 1,... 4, at four values of the 
light-matter coupling strentgh. The lower panels show that at large magnetic fields (B = 10T) both the m = 3 and m = 4 LL 
the pseudospin polarization still remains degenerate. Hence higher order perturbation terms should be necessary in order to 
lift this degeneracy. 


of experimentally realizing such a superposition state. A 
paradigmatic case of such interesting superposition states 
is given by the coherent state which are minimal uncer¬ 
tainty wave packets relevant for studying the classical 
states of the radiation field in the sense of being a classical 
counterpart of the quantum harmonic oscillator. Since, 
as in the quantum harmonic oscillator, the Landau levels 
in graphene are eigenstates of the number operator, one 
should expect that these coherent state superpositions 
would be interesting. Indeed, there has been recently 
some proposals to analyze the dynamics of coherent elec¬ 
tronic states in graphene nanomechanical resonators^. 
By taking advantage of the intrinsic non linear nature 
of flexural modes in graphene, the authors of reference^! 
show that cat-like^ states can be generated. We follow a 
different physical approach. Instead of resorting to non- 
linearities of flexural modes we invoke the light-matter 
coupling as a mechanism for studying the evolution of 
an initially prepared coherent superposition state built 
from the Landau level eigenstates described in section II 
and evaluate the induced pseudospin polarization effects 


in order to contrast to the results shown in FIG. [2] 
Formally speaking, the coherent state jcu) is defined by 
means of the eigenvalue equation 

A\a) = a\a ), (40) 

where A = a 1. Using the expansion 

I <*> =c 0 |<A)) (41) 

sn 

with X)sn representing a summation for all n ^ 0. As¬ 
suming, without loss of generality, the symmetric sce¬ 
nario c~ s = c®, one finds that the coherent state is given 
as 

|a> = e-^ (Vo) + 4E^=bn)) > (42) 

which can be shown to be normalized. If we now evalu¬ 
ate the mean value of the hermitian operator A^A in the 
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FIG. 3. Time dependence of the pseudospin polarization as given in equation (45). The panels show the dynamical behavior at 
four characteristic values |a|. For the numerical evaluation of the series we have truncated at n = N = 100. The upper panels 
show the coherence among the lowest LL whereas the lower panel show localization effects approaching the classical behavior, 
which corresponds to large values of |a|. See the discussion in the main text. 


coherent state and use the definition given in equation 
mop , we get (a|AtA|a) = \a\ 2 . In addition, if we use 
the coherent state to evaluate the average of the number 
operator N a defined in equation (TIBI) , it is not difficult 
to show that we also get (a\N a \a) = \a\ 2 . Therefore, the 
coherent state parameter amplitude |a| is a measure of 
the mean number of Landau levels that are excited and 
corresponds to the mean photon number in the context of 
quantum optics. Thus, we will use it as a control param¬ 
eter to discuss the properties of the physical quantities 
as follows. 

In this case, we find for the pseudospin polarization 


{a\a z {t)\a) 



n 


(hu c ) 2 cos 2 e n t/h+ n£ 2 


( huj c ) 2 + n £ 2 


(43) 


This is plotted in FIG. 3 for £ = 0.25 tko and four rep¬ 
resentative values of |a|. In this figure we notice that 
for small values of \a\ the dynamics of the pseudospin 
polarization resembles the pattern for Rabi oscillations 
since the main contributions would arise for the interfer¬ 
ence among the zero and first LL. Yet, no polarization 


inversion can be achieved within this regime. 

However, once |a| = 1 the contribution from other 
LL states becomes increasingly important to the inter¬ 
ference pattern and the former Rabi oscillations become 
distorted. Moreover, at this value of |a|, one can achieve 
the polarization inversion for large enough values of time 
in the long-term evolution. In the lower panels of FIG. 
4 we find that for larger values (|a| = 4 and |a| = 8 ), 
we get a beating pattern showing a dynamical localiza¬ 
tion effect that is directly related to a collective behavior 
of the driven charge carries in graphene. We would like 
to remark that the problem of the population inversion 
has already been studied long time ago by Eberly et &&. 
Their model corresponds to a two-level system coupled to 
a single mode radiation field (Jaynes Cummings Hamilto¬ 
nian). In our approach, this in turn is given by the static 
graphene Hamiltonian for Landau levels written in equa¬ 
tion (2). A larger value of the coherent state parameter 
| a | would imply larger mean Landau level occupation, 
and thus the second term in the numerator of the pseu¬ 
dospin polarization given in expression (43) would have 
a larger influence in the pseudospin inversion since larger 
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FIG. 4. Time dependence of the autocorrelation function as given in equation (48) truncated at n = N = 100. For large |a| = 4 
and |a| = 8, it is the cyclotron frequency oj c instead of the driving frequency c o which determines the time scale for the revival 
times. Comparing to FIG. 3 it is apparent that information on the pseudospin dynamics can be indirectly inferred from the 
autocorrelation function dynamics. See discussion in the main text. 


n Landau levels would be occupied which will contribute 
a higher weight in the pseudospin polarization (see nu¬ 
merator in equation (20). Moreover, since the effective 
light-matter coupling strength £ also affects the phase, as 
given by the cosine term in equation (20) of the polariza¬ 
tion, this interesting interplay forbids the inversion which 
is also a signature of the localization and beating effects 
shown on panels 3(c) and 3(d) and our results extend 
those found in the context of reference^. 

This physical picture for the coherent state dynamics can 
be complemented by studying the autocorrelation func¬ 
tion C a (t) = (a'|\f , (f)), which is found to be given as 

C a (t ) = e - ^ ^ j cos e n t/h + i— sine n t/h j . 

n n - \ tn J 

(44) 

The autocorrelation function provides additional physi¬ 
cal information of the system since its Fourier transform 
is related to the local density of states^. Its time evo¬ 
lution is plotted in FIG. 4 choosing again an effective 
coupling strength value of £ = 0.25 Hu> and the same val¬ 


ues of a as in FIG. 3 showing the dynamical behavior 
of the pseudospin polarization. Comparing FIG.3 and 
FIG. 4 we see that partial revivals for the autocorrela¬ 
tion function are correlated to the beating or localization 
behavior of the pseudospin polarization. Therefore, one 
could indirectly gain information on the pseudospin dy¬ 
namics by measuring the time revivals^ which means 
detecting those times for which the wave packet recon¬ 
structs itself. In order to show explicitly the role of the 
driving field we have plotted in FIG. 5 the static auto¬ 
correlation function. Comparing the lower right panels 
in FIG. 4 and FIG. 5, we find that the radiation field 
induces additional revivals in accordance to the beating 
pattern in the pseudospin oscillation, as seen in FIG.3. 


IV. DISCUSSION 

We have shown that the radiadion field leads to a 
quasienergy spectrum with a level dependent gap, ex¬ 
cept for the to = 0 LL which, to leading order, remains 
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FIG. 5. Time dependence of the autocorrelation function as given in equation (48) truncated at n = N = 100 in absence 
of radiation field £ = 0. In the upper panels we see that at small values of |«| = 0.25 and |a| = 1, the static and driven 
autocorrelation function are qualitatively similar. Yet, the photoinduced quantum revivals at larger times are not seen in the 
two lower panels as compared to FIG.4 


insensitive to the radiation field effects. In addition, we 
have found that a finite out of plane pseudospin polariza¬ 
tion value can arise for initial states that possess either 
a finite or vanishing initial value of er z . As we saw at 
the beginning of section dm in presence of the radiation 
field, the total out of plane angular momentum compo¬ 
nent j z of the electrons is no longer a constant of motion. 
Therefore, the pseudospin oscillations are due to the an¬ 
gular momentum exchange among the driving field and 
the charge carriers in single layer graphene. For the co¬ 
herent state a beating pattern emerges in the pseudospin 
polarization. In spite of the fact that the interaction 
is dictated by a periodic Hamiltonian we find that the 
dominant or characteristic time scale for this collective 
behavior is the cyclotron frequency. We could expect 
this behavior of the wave packet at large values of |a| as 
a measure of the classical behavior associated to the cy¬ 
clotron problem but now realized with Dirac fermions in 
the LL quantized regime as was discussed by one of the 
authors in reference^. We must remark that in order to 
be able to detect the reported effects an ensemble of co¬ 
herent states should be prepared each time, because the 


measurement process destroys the state. It should also 
be remembered that the coherent state parameter |a| is 
determined by choosing appropiate values for the mean 
values of position and momentum in accordance to the 
following prescription 

{a\x\a) = (45) 

{a\p x \a) = ^-S(a). (46) 

Another point to be highlighted is that for standard two 
dimensional electron gases the coherent state built from 
the Landau levels would remain coherent, i.e., it will 
evolve in time in such a manner that it would just os¬ 
cillate in time around the prescribed mean values given 
in eqs. (14511 and (l46l) . More precisely, its Wigner function 
representation in phase space will oscillate in time with¬ 
out deformation. This is known to be a consequence of 
the fact that the dependence with the n quantum num¬ 
ber is linear. Yet, in graphene we have a y/n and this 
in turn prevents the coherent state to evolve in such a 
coherent manner. 
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Therefore, the spreading of the wave packet and the 
corresponding appearance of the additional revival times 
at shorter time values is a direct consequence of the driv¬ 
ing field that even at low coupling can induce interesting 
dynamical behavior as shown in FIG. 3 and FIG. 4. We 
supported this last statement by evaluating the static 
(£ = 0) autocorrelation function as shown in FIG. 5 and 
found that the second packet reconstruction (see lower 
panel for |ck| = 8 in FIG. 4) around tu c = 24 is absent in 
the static regime. 

We would also like to comment that graphene subject 
to electromagnetic radiation without quantizing mag¬ 
netic field has been discussed in several papers^. It has 
been shown that for zero momentum (k=0), the dynam¬ 
ical equation is exactly solvable and leads to a photoin- 
duced (zero quantizing magnetic field) mass term that 
has opposite signs at the two valleys. Yet, at finite mo¬ 
mentum the dynamics is no longer exactly solvable and 
one has to resort to numerical analysis by means of an 
infinite expansion in Fourier modes. This is why we con¬ 
sider our results to be a valuable tool in analytically de¬ 
scribing the photoinduced pseudospin effects in the Lan¬ 
dau level structure of monolayer graphene. 

We now briefly mention that our values of quantizing 
magnetic fields are within experimentally accesible or¬ 
ders. For instance, in the pioneering paper by Novoselov 
et al. in reference^, the authors used a value of quantiz¬ 
ing magnetic field B = 14T to study the quantum Hall 
effect in graphene. Moreover, in a following classical pa¬ 
per, their results were extended up to values of B = 29T 
and even B = 45T for the quantizing magnetic fields. 
Here it is shown that at these values of magnetic field, 
the quantum Hall effect in graphene could be observed 
at room temperature that constitutes a breakthrough in 
the physics of quantum Hall phases. 


V. CONCLUSIONS 

We have analyzed the dynamical modulation of 
physical quantities for Dirac fermions within the Landau 
level quantized regime of single layer graphene subject 
to intense circularly polarized Terahertz radiation. To 
the best of our knowledge, this is the first analysis of 
the photoinduced manipulation of the LL structure in 
single-layer graphene subject to a continuously applied 
intense laser field instead of a pulsed one. By means 
of a perturbative analytical treatment we found very 
interesting physical features such as a non trivial level 
dependent dynamically induced gap structure. Due to 
the angular momentum exchange among the radiation 
field and the charge carriers, it also leads to modulation 
of the oscillations in the dynamics of the out of plane 
pseudospin polarization, even for superposition states 
with an initially vanishing pseudospin value. We also 
found that localization effects in the time evolution of 
the pseudospin polarization can be keep track of by 
measuring the revival times of a wave packet initially 


prepared as a coherent state. The reported photoinduced 
gap modulation and pseudospin oscillations could be 
detected through the reemitted dipolar radiation from 
the oscillating charge carriers as it was proposed in 
reference^. 

We would also like to mention that for values of A 1, 
one enters the so called ultrastrong coupling limit in the 
context of cavity quantum electrodynamics experiments. 
Within this parameter regime, it has been experimen¬ 
tally shown in reference^! that explicit anti-crossings ap¬ 
pear in the energy spectrum. This in turn corresponds 
to a breakdown of the Jaynes cummings approximation 
that describes resonant processes in the two-level prob¬ 
lem coupled to a single mode radiation field. Our exact 
Floquet Hamiltonian given in equation (20) allows us to 
explore this ultrastrong coupling regime beyond the per¬ 
turbative results presented and could be the subject of 
future work where one would expect non linear effects to 
be relevant in the graphene physics. For instance, as we 
already argued in the discussion section, one could ex¬ 
pect that the accidental degeneracies seen in the panels 
of FIG. 2, at large values of the quantizing magnetic field, 
would be lifted by a stronger coupling to the radiation 
field. 

Finally, we would like to mention that although we per¬ 
form our analysis for neutral undoped graphene. Yet, 
in the graphene literature the role of doping has been 
given special attention since it might either modify the 
effective Fermi velocity of Dirac fermion s 28 ! 29 . The first 
workSS shows that the energy spectrum near the charge 
neutrality point is non linear and no gap is found at en¬ 
ergies even as close to the Dirac point as O.lmeV. On 
the other hand, the authors of paper— show that that 
electron interactions leave the graphene energy disper¬ 
sion linear as a function of excitation energy even for 
energies within ±200meV of the Fermi energy. 

Yet, it has recently been shown that doping effects might 
also lead to a bandgap opening in the graphene spec¬ 
trum (see^H and references therein). In particular, the 
authors of the worb22 have discussed the tunability of 
the bandgap energy in single-layer graphene due to man¬ 
ganese oxide nanoparticles by means of an electrochem¬ 
ical method^. This work reports a maximum value for 
the induced energy bandgap of 0.256 eV. Although they 
do not discuss the quantized LL regime, we could in¬ 
corporate these doping effects in our model by including 
a phenomenological diagonal term in the Landau level 
Hamiltonian given in equations (2) and (3). This mass- 
term would in turn be proportional to a z and the valley 
index rj, as required by time reversal invariance. Thus, 
one would expect an interesting interplay between the 
photoinduced gap openings and the static bandgap de¬ 
termined by doping effects. Yet, we considered undoped 
graphene in order to highlight the photoinduced LL de¬ 
pendent bandgap discussed above. 

Concerning the actual experimental observation of our 
predicted results, we consider that these could motive 
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the exploration of larger intensities for the radiation field 
as given in the experimental setup in reference^ when 
their model enters the Landau level regime. In this sense, 
we consider that our work would contribute to explore a 
new physical scenario within realistic parameter values 
to discuss the pseudospin physics in graphene, also 
taking into account the role of Zitterbewegung as it 
is discussed in reference^, extending their results to 
initially prepared coherent states when the radiation 
field is a monochromatic continuous laser field instead 
of a pulsed laser. 
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VI. APPENDIX 

A. Perturbative calculation of effective 
Hamiltonian for single layer graphene 

In order to get the effective Hamiltonian for single layer 
graphene we need to evaluate the following expression 

H = e x/ 2 I -H F e~ x/2I ~, (47) 


with the antihermitian operator 

I- = aV_ — aa + . (48) 

Using the Baker Campbell Hausdorff formula we have 

H = H F + ±[I-,H F ] + ±(?) [I_,[I_,H f ]} + ... (49) 

The first commutator is worked out explicitly 

[/_, H F ) = [a^cr- — dc+, H F ], (50) 

Since H F = w c / + — ojN a + £ a x and [I-,N a ] = 0, with 
1 + = cAa- + ao+, we only need to evaluate two commu¬ 
tators. The first one gives 

[/_,/+] = [a^<7_ — acr+,dV__ + acx+] 

= — 2[dcr+, a)<j-] 

= —2(a[cr+, aV_] + [a, dV_]cr + ) 

= — 2(aal[cr+, <r_] + [a, a^]cr_cr+) 

= —2(a a)(T z + (1 — cr z )/2) 

= —2(d^acr z + (1 + CTz)/2) 

= -2(a) a + (1 + <r z )/2)cr z 

= -2N a a z , (51) 

whereas the second one follows as 

[J_, CT X ] = [dV. - dcr+, <7_ + cr+] 

= — d[tr + , er_] + [cr_, cr+] 

= — (a + a^)cr z . 

Upon substitution of these first-order corrections and in¬ 
troduction of the shifted harmonic oscillator operators 
we get the effective Hamiltonian given in equation (1231) . 
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